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Abstract 

Upon initial excitation of a few normal modes the energy distribu- 
tion among all modes of a nonlinear atomic chain (the Fermi-Pasta- 
Ulam model) exhibits exponential localization on large time scales. At 
the same time resonant anomalies (peaks) are observed in its weakly 
excited tail for long times preceding equipartition. We observe a sim- 
ilar resonant tail structure also for exact time-periodic Lyapunov or- 
bits, coined q-breathers due to their exponential localization in modal 
space. We give a simple explanation for this structure in terms of 
superharmonic resonances. The resonance analysis agrees very well 
with numerical results and has predictive power. We extend a previ- 
ously developed perturbation method, based essentially on a Poincare- 
Lindstcdt scheme, in order to account for these resonances, and in 
order to treat more general model cases, including truncated Toda po- 
tentials. Our results give qualitative and semiquantitative account for 
the superharmonic resonances of q-breathers and natural packets. 
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The Fermi-Pasta-Ulam model (FPU) is a one dimensional 
nonlinear lattice whose dynamics, for small energy densities £ 
and initially localized (in Fourier space) excitations, is charac- 
terized by a long time freezing of the harmonic energy among 
a subset of normal modes (natural packets). This metastable 
regime persists until, after times growing exponentially in j, 
equipartition is reached. The intermediate state is character- 
ized by an exponential localization of the harmonic energies 
around the initially excited modes and the appearance of res- 
onant deviations in the tail of the energy distribution. During 
the transition from metastability to stochasticity, redistribu- 
tion of mode energies happens through the tail resonances. 

The system admits also exact periodic orbits - q-breathers 
(QB) - which are obtained by continuation of the normal 
modes of the linear system. QBs also show both exponen- 
tial localization of mode energy distribution and resonant tail 
structures. The QB localization holds below an energy den- 
sity threshold £ * which depends on the nonlinear parameter. 
The tail resonances are due to small denominators which ap- 
pear naturally in a Poincare-Lindstedt perturbation scheme. 

We observe that the dynamics of the nonlinear FPU model 
generated by initially exciting one normal mode (FPU tra- 
jectory) exhibits the main elements of the nearby g-breather. 
Thus the presented analytical approach, besides describing 
the g-breather, also explains some relevant features of the 
corresponding FPU trajectory. Remarkably, this correspon- 
dence may hold on very long times, which requires further 
rigorous explanations. 



1 Introduction 

1.1 Brief history 

The computational experiments which Fermi, Pasta and Ulam per- 
formed around 1950 in the Los Alamos laboratory are nowadays con- 
sidered as an important step in the progress of understanding nonlinear 
dynamics. The main aim of their numerical investigation was to ob- 
serve, as they expected, a "relatively fast" and "equally distributed" 
sharing of the averaged harmonic energy of an anharmonic atomic 
chain among all the normal modes, even when starting with initial 
conditions far from some equipartition regime. When chosing initial 
data localized on just one low frequency mode they however did not ob- 
serve any actual equipartition, but a localization of the energy among 
some initial modes in the nearby frequency sector (also coined in the 
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literature and hereafter natural packets), including some additional 
almost complete recurrence to the initially excited mode. 

Among the several explanations that followed, we recall the one 
given in 1965 by Izrailev and Chirikov |23l I24| . They suggested the 
idea that the KAM theorem was the right approach in order to explain 
the regular dynamical behaviour observed: if the energy given to the 
system is below a threshold E* (N) , then the dynamics is influenced 
by the persistence of most of the nonresonant invariant tori. Since 
lmi;v->oo£* = 0, where £* = E j s the energy density the FPU 
phenomenon should not be relevant in the thermodynamic limit at 
any finite temperature or energy density. In 1970, instead, Bocchieri, 
Scotti, Bearzi and Loinger [S] conjectured the existence of a specific en- 
ergy threshold £* (which remains positive in the thermodynamic limit) 
below which equipartition is not reachable. Recent numerical results 
OEEZj indicate that below that threshold equipartition is still reached 
for infinite times, yet the time scale of the route to equipartition may 
diverge exponentially with lowering the energy density. 

Several computational experiments have been performed (e.g. |111 
E3 E3 |3 H ED), suggesting possible answers about the time needed 
to approach equipartition, its scaling with the energy density and di- 
mension of the system and the features of the dynamics during this 
transition. Only few analytical studies touch the case of large systems: 
here we mention some relevant results related to the application of the 
KAM Theorem [311 133j. to some special (periodic) solutions |291 132j 
and to a normal form construction via a PDE 1 . 

Recently Flach, Ivanchenko and Kanakov ^H] have focussed on the 
main FPU observation that the initially excited normal mode shares 
its energy for long times only with a few other modes from a frequency 
neighbourhood in modal space. They have identified this long lasting 
regime as a dynamical localization effect and applied the methods de- 
veloped for discrete breathers in FPU chains ^S] to the dynamics of 
normal modes. The result is that q-breathers (QB) - time-periodic and 
modal-space-localized orbits - persist in the FPU model. The dynam- 
ics generated by one initially excited mode evolves close to the related 
q-breathers for very long times. Thus many features of the short- and 
medium-time evolution of natural packets are encoded in the profile 
of a QB and the phase space flow of small fluctuations around a given 
QB. 

1.2 The FPU model 

Fermi, Pasta and Ulam considered a non-linear one-dimensional lattice, 
with fixed boundary conditions 

[ H(x,y) = | Y,f=i v) + E^Lo V{x j+1 ~ x i)> (i) 
\x = x N+1 = 0, 

where V(r) = \r 2 + jr 3 + jr 4 is the coupling potential, (Xj,Uj) are 
respectively the displacement from the equilibrium of the j-particle and 
its kinetic momentum, while the parameters a, (3 may be considered 
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as real and positive. The case of linear equations of motion, given by 
the quadratic Hamiltonian 

1 N 

H (x,y) = - [Vj + (xj+i - xj) 2 ] , 
3=0 

is integrable and corresponds to a system of TV uncoupled normal mode 
oscillators. We can show this via the change of coordinates 




which diagonalizes the Hamiltonian Hq 

N 

2 



i N 



where (Qj,pj) are the normal mode coordinates and momenta, and i 
are the normal mode frequencies given by 

• ( j n 
Qjj = 2 sm 

3 \2N + 2 

The effect of the non-linear part in the model |QJ is to introduce an 
interaction between the different normal mode oscillators. In this way 
they exchange energy and the harmonic actions are in general no longer 
preserved. Even if the considered change to normal modes appears 
hard to develop for the cubic and quartic part of the Hamiltonian, 
there is a mathematical way to exhibit which are the selection rules 
that couple the different oscillators (see also QOl)- Let us consider the 
so called a-model, where = 0. The cubic part of the potential can 
be rewritten as 

JV 

3 ( x 3+ l ~ = 3 ~7== T 2 ^ WiUjUhQiQjQhBijh, 

j=0 l<i,j,h<N 

where we set B ijh = X)<r 1>2 =±i ^i+Tij+T 2 h, with 

,i±j±h = 

&i± 3 ±h={-l ,i±j±h = 2(N + l) . (2) 
, otherwise 
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Figure 1: (a): Natural packet evolution initially exciting mode qo = 1 
The energies of normal modes are plotted versus q. Circles - 10 4 , squares • 
10 5 , rhombs - 10 6 . Dashed line - g-breather from (b) for comparison, (b) 
The energies of normal modes versus q for a g-breather with qo = 1. (c) 
sequence rjj . (d): sequence Xj ■ Parameters in all cases: N = 31, a = 
0.33, £ = 0.01. 



The equations of motion for a normal mode oscillator then read 



7T UJjUliUJhQiQhBijh. 
1/ Ki,h<N 



(3) 



1.3 Natural packets 

One of the findings which can be related to Bocchieri's conjecture is 
the so called 'metastability' of natural packets, which has been nu- 
merically investigated by A.Giorgilli and others EI El EZj • What 
has been observed in [3] is that, for energy density small enough, an 
initial excitation of the first mode is spread among nearby modes up 
to a natural frequency which depends only on the energy density, ac- 
cording to the power law u> ~ f 1 / 4 . This localization of the harmonic 
energy on a subset of modes, which can coexist with chaotic dynam- 
ics inside the subset, persists for times which can grow exponentially 
with the inverse of the energy density (see 0]), and only after this long 
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Figure 2: (a): Natural packet evolution initially exciting mode qo = 2 
(circles). The energies of normal modes are plotted versus q at t = 10 5 . The 
dashed line shows the energy distribution for a corresponding g-breather 
with qo = 2. (b): Same as (a), but for qo = 3. (c): sequence ry~ 1 (2,47). (d): 
sequence ry^ 1 (3, 47). Parameters in all cases: N = 47, a = 0.1, £ = 0.0025. 
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transient the system reaches a proper equipartition regime. The time- 
dependence of the size of these natural packets of normal modes show 
two time scales n,2 which separate strongly when the energy density 
is below a certain threshold. The first time scale t\ characterizes the 
time needed for the formation of a natural packet of nearly fixed (and 
quite localized) shape in normal mode space. The second time scale t%, 
which becomes much larger compared to the first one when the energy 
density is lowered, characterizes the time needed to finally delocalize 
the packet over all modes and to reach equipartition. Both time scales 
Ti « T2 at some energy density threshold. Thus at and above the 
energy density threshold no metastable or intermediate normal mode 
localization occurs, and the natural packet spreads more or less con- 
tinuously in time over all modes. Further numerical simulations 0] 
used the spectral entropy indicator (which gives the effective fraction 
of modes involved in the dynamics) and showed that the second time 
scale nicely fits with an exponential law. 

An example of the evolution of such a packet for energy densities 
below the threshold (when n < 10 4 , t-i > 10 9 ) was observed in [2] 
(Fig. 3 therein). Following that work, we performed a similar simula- 
tion, see Fignja). The normal mode qo = 1 has been initially excited, 
and the harmonic energy distribution is shown for three different times 
t = 10 4 , 10 5 , 10 6 . Since harmonic mode energies may fluctuate, we plot 
the time-averaged quantities, with a time-averaging window of 10 4 . 
The logarithmic scale reveals that for t — 10 4 the overall strong expo- 
nential distribution, which is locked over many decades in time, shows 
pronounced resonant peaks in the far tails (note however that even in 
the peak maxima the mode energies are orders of magnitude smaller 
than the energies of the main modes excited). At t = 10 5 a gradual 
change is observed - the resonant peaks in the tail slowly grow and 
drag the tail upwards. At t = 10 6 , approaching r 2 , this process con- 
tinues and the slow but steady spreading of the packet over all modes 
can be already anticipated. At even larger times, the harmonic energy 
becomes equally distributed over all the normal modes (not displayed 
here). 

In Fig. H[a,b) instead, the initial excitation is localized on the sec- 
ond and third mode, respectively. The profiles arc shown for time 
T\ < t = 10 5 < T2- They are again characterized by the existence 
of peaks in the tails of the distributions. The persistence of these 
anomalies when changing the initial mode number excited, suggests 
that there is a robust and generic pathway of equipartition which goes 
via excitation of these resonant tail peaks. 

Let us compare these results to the profiles of corresponding q- 
breathers ^H]- We compute periodic orbits, which correspond to only 
one normal mode with qo being excited in the limit of vanishing non- 
linearity. In Fig^Jb) the profile of a g-breather with qo = 1 is shown. 
The harmonic normal mode energies are averaged over one period of 
the g-breather. We observe exponential localization, and in the tail two 
distinct resonant peaks. For comparison we plot the same result as a 
dashed line in Fig^a). The exponential profiles of the q- breather and 
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the corresponding natural packet are strikingly similar in the whole q- 
space. That supports previous results that a natural packet for t <C 
evolves very close to a corresponding g-breather orbit. The tail reso- 
nance locations though differ for both trajectories. That signals that 
the nature of the tail resonances can be understood already by study- 
ing a g-breather, but their location is perhaps very sensitive to small 
perturbations. Nevertheless some common features are also observed 
- the tail resonances start to appear in the second part of the q-space. 
The dashed lines in Fig|2Ja,b) are g-breather profiles for corresponding 
smaller energies. We note that here not only the exponential decay fits 
well with the natural packets, but now also the location and even the 
strength of the tail resonances agrees well with the one of the natural 
packets. We conclude that the correspondence between natural packets 
and g-breathers improves with lowering the energy also with respect 
to the sensitive tail resonance structures. 

Further similarity between natural packets and g-breathers follows 
from [22j where it was shown, that the localization length of a QB de- 
pends on the expansion parameters of the potentials, the energy density 
£ and the normalized wave number ko = qo/N, but not anymore on 
the system size N. In fact QBs can be scaled to and obtained in in- 
finitely large lattices [2HI • Using estimates from perturbation theory we 
can conclude that delocalization of QBs will happen when the energy 
density exceeds some critical value (with the other relevant parameters 
being fixed). It can be expected then that an initially localized excita- 
tion will quickly spread over all modes as well. This is in contrast to 
the case when the energy density is below the QB threshold. Then the 
QB is exponentially localized, and an initially localized excitation is 
known to spread a bit and take a form similar to a QB for long times, 
with subsequent departure for much longer times. That gives further 
evidence that the exact and simple periodic orbits - q-breathers - de- 
scribe the metastable features of a natural packet spreading in modal 
space. It is then also natural to study whether the q-breather shows 
similar resonant tail peak features as the natural packet which evolves 
close to the QB. For example, in Fig.^b) the exact QB corresponding 
to go = 1 shows an exponential decay which persists also in the tail, 
while the tail's harmonic distribution of the close orbit (Fig. [Ja)) is 
more irregular. Moreover, the QB also exhibits, even if in different 
sites, a couple of resonant peaks. 

In the following we will first (Section 2.) formulate a simple super- 
harmonic resonance analysis and show that it accounts quantitatively 
for the observed tail peak positions. Then (Section 3.) we will perform 
a modified perturbation approach to the q-breather profile and provide 
a semiquantitative answer for the tail profiles. 
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2 Tail resonances 



2.1 g-breathers 

As we noticed concluding the previous section, a fruitful way to inves- 
tigate the dynamics for very localized initial excitations is to study the 
existence and the related properties of nearby invariant manifolds, we 
thus focus first on the Lyapunov orbits, called g-breathers due to their 
localization in the normal mode space. These orbits are the natural 
continuation, for small values of a coupling parameter, of the uncou- 
pled oscillatory solutions which appear in the linear problem [lfcij . The 
existence of these families for every linear frequency u)j of goes 
through classical continuation theorems, such as the Lyapunov Theo- 
rem or the Poincare Theorem, and essentially relies on the property 
of rational independence of any couple of frequencies composing the 
spectrum. Indeed old (0) and more recent (pU) results on sums of 
roots of Unity have proved that any linear combination of two different 
frequencies with rational coefficients is always different from zero 

kjujj + km ^ 0, Vfcj,fc;GQ, luj ^ wj. (4) 

In order to apply the above mentioned Theorems, it is enough to use 
a weaker result which ensures that Q holds when ki = 1 and u>j is 
the frequency of the linear periodic orbit we want to continue. The 
new periodic orbits are expected to be continuous deformations of the 
original ones upon changing the nonlinearity parameter. Thus, if we 
initially excite only one go _m ode, we can imagine the dynamics be- 
ing strongly influenced by the corresponding go-breather: in fact, the 
initial condition is close to this special orbit in phase space. 

The idea to interpret the original results obtained by Fermi, Pasta 
and Ulam in terms of deformation of periodic Lyapunov orbits has 
been recently used in [lfi| . In this paper the authors perform both a 
numerical and analytical investigation of the a and /3 models. Their 
study gives an exponential localization of the energy among the modes 
involved in the periodic orbit; these modes are selected according to 
the definition of the coupling coefficients Bijh (see eq (0)). Further- 
more, the linear stability analysis for the fi model remarkably reveals 
that this exponential localization, which gets weaker when increasing 
the coupling parameter /3, is still persistent when the orbit becomes 
unstable. From the analytical point of view, the authors provide an 
application of the century old proof of Lyapunov to the existence of 
q-breathers, and a perturbative Poincare-Lindstedt expansion which 
explains the localization for the modes kqo belonging to the low fre- 
quency part of the spectrum, id est when both qo <C N and kqo <C N. 
This perturbation scheme also shows how, at every step fc, only one 
new mode, precisely the mode kqo, is involved; then, by using a cubic 
expansion of the spectrum to cope with the small divisors, the dynam- 
ics of this mode is approximated as an oscillation whose frequency is 
kuj qo . 

In order to better understand the idea that moved our research, we 
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recall here just the main steps of their deduction, and will extend this 
approach in the next section. The starting point consists in expanding 
the periodic solution with respect to a perturbation parameter, here 
chosen as a = 

Q q {t)=Y t °*Qf® ( 5 ) 

i>o 

in |JSJ. Sorting the terms with respect to powers of a, we get an infinite 
sequence of differential equations, one for each component Q^ q \ j > 0; 
the zero order solution is simply the cosinusoidal oscillation of the qo 
mode (or the linear approximation) 

<2g 0) = <*9,*Ao COS (uj qo t). 

For any n > 2, the scheme then follows these rules: 

• at every step n — 1, only one new mode Qlq ^(t) is involved in 
the motion. This is related to the coupling coefficients B nqQt i_ m . 

• up to the step n — 1, every solution Q q h \ < h < n — 2,qo < 
q < (n — l)qo is a finite sum of harmonics having frequencies 
muj qo , m < h+1. Among these solutions, the one Q% ^{t), j = 
2, . . . , n — 1 are approximated with the largest amplitude oscilla- 
tion. Therefore the new forced oscillator becomes 

j2 n 

^Qnqo + ^IgoQnqo = E C m COS (mcu qo t), 

m=0 

with a suitably defined sequence C m ; the initial conditions are 
chosen such that the solution reads 

n „ 

= E U 2 .l,^ COS Kt). (6) 
m =0 "90 Ul W qo 

• exploiting the cubic expansion at low frequencies LOk ~ k5 — 
^k 3 S 3 , S = jfrj, it is possible both to realize that the small- 
est divisor is o;„ go — n 2 ui 2 o and to give an approximated value of 
this quantity as 



< -n 2 <«-^V-lK, (7-. 



Hence, the solution can still be reduced to a single oscillator 
(n-i) 12C n i +\ ( ~\ /o f n 1o\ 2 



if the error (^f) 2 is small. 
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Figure 3: Parameters: N = 47, a = 0.1, £ = 0.0025. (a): Energy of normal 
modes versus mode number of a QB for go = 2. (b): Similar to (a), but for 
a QB with qo = 3. (c): sequence xj 1 (2,47). (d): sequence Xj 1 (3, 47). 



2.2 Higher order nonlinear resonances 

The perturbation scheme used in is based on the cubic approx- 
imation of the dependence of Wj on j. This approximation holds in 
the low frequency domain and allows to perform a detailed analysis of 
the small divisors. In this way one limits the study to the first third 
of the spectrum from below, were this approximation holds. However, 
the numerical computation of these orbits clearly shows an exponential 
decay which continues into the upper two thirds of the normal mode 
frequency spectrum. However, some of these orbits also show resonant 
deviations in the tails. We conjecture that these resonances are due 
to almost resonances which are not captured by the abovementioned 
cubic approximation of the dependence of normal mode frequencies on 
mode numbers. Thus we investigate the sequence of small divisors 

Vn,m = k 2 i«o _ m2uj q 1 m = n,n-2,n-4,... (8) 

for every frequency U! nqo belonging to a mode involved in the (/-breather 
considered. The special choice of only odd or even multiples is related 
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to the perturbation theory and will become clear later. Then, for each 
n, we can also consider the sequence of the smallest divisors 

Vn = min {rj nm }. (9) 

m—n,n—2,... 

For each qo, the above definitions refers to u) qo , which is the first order 
approximation of the true frequency tt qo of the go-breather considered. 
Equations |(HJ| and @ are a consequence of applying the Poincare- 
Lindstedt method, since these are the small denominators appearing 
in the Fourier expansion of the leading solution JJjJ. However the ap- 
proximating fi go by oj qo is less accurate the larger we take the product 
of perturbation parameters a£. Thus intuition suggests to consider 
also sequences obtained by replacing u qo with the true QB frequency 
Sl qo {a,£) 

X n , m = \ul qo ~ m 2 fl 2 \, X n= min {x„. m }, (10) 

' y H> m— n,n— 2,... ' 

with the same choice of m. These are the small denominators appear- 
ing in a pure Fourier expansion of the solutions Q n q a by substituting 

Qnqa = EkezQnq ,ke lknmt HI © . 

In Fignjc,d) we plot the sequences 77" 1 and \n 1 f° r Qo = 1) N = 31 
and £ = 0.01, and in Fig. 01c,d) and Fig|3fc,d) respectively the se- 
quences r)- 1 and X n l for q a = 2, 3, N+l = 48, a = 0.1 and £ = 0.0025. 
Note that the corresponding true QB frequency fl qo has been used to 
compute X- I n au cases the first part of the sequence is influenced by 
the cubic approximation. The curves follow the law 3 \ g(kS) an( ^ ^ ne 
largest term corresponds to the smallest allowed value of q. Far from 
the low frequency sector, we can observe a different and irregular be- 
haviour. It happens to some of the frequencies to be strongly almost 
resonant with multiples of the QB frequency uj qo (or f2 9o ). Let us es- 
timate the lower bound k c qo for the existence of such strong almost 
resonances. It is enough to find the first k c for which 

w fcc90 - (fe c - l)uj qo = 

holds. Expanding uj q for small q to third order we find k c sa (^ff^)i , or 

in size-independent quantities k c ~ (^p?) 3 • This value agrees well with 
the location of the first local minimum in the mentioned figures. The 
subsequent resonances suggest that the QB profile will show resonant 
deviations in its tail, exactly at the predicted locations. The sequence 
Xj for very small nonlinearities can be well approximated by In- 
creasing the perturbation, new small denominators may appear and 
some of the old one can become less relevant due to the QB frequency 
change. 

In Fig. ^d) the x sequence corresponds to the QB in Fig. ^b). 
The resonant local minima for Xj are J = 16 and j — 21, which per- 
fectly agree with the two tail peaks of the corresponding QB. Note 
that the ?y sequence in Fig. QJc) shows almost the same resonance loca- 
tions, but the relative strength is different. Also both sequences show 
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additional resonance locations, which are not observed in the QB tail. 
Instead almost all of them are observed in the energy distribution of 
the corresponding natural packet in Fig. HJa) . 

The distributions in Fig. 0{a,b) refer to two natural packets local- 
ized respectively on qq = 2 and go = 3, while in Fig. EIa,b) we show 
the corresponding exact go-breathers, the profiles, including the (much 
weaker) tail resonances, coincide remarkably. That is due to the weak 
nonlinearity. The corresponding 7/ sequencies in Fig. |2Ic,d) and x se ~ 
quencies Fig. |3Jc,d) also show very good agreement with the really 
observed locations of tail resonances. 

A conjecture which explains the existence of the observed peaks 
in the tails can be now formulated. The amplitudes of some modes 
are amplified by very small denominators generated via strong almost 
exact resonances of the type ui^ qo — m 2 rt qg . The -q and \ sequencies 
allow to predict possible locations of these supcrharmonic resonances, 
however not always the predicted locations are indeed occupied by 
observed resonances. 

Let us now formulate a set of QB properties which have to be stud- 
ied. First we recall again that the energy distribution of the computed 
g-breathers for go <C N shows that their exponential localization is 
kept throughout the whole g-space and is not restricted to the low q 
mode sector, as anticipated by the approximations used in perturba- 
tion theory. Hence, it is an important task to extend the perturbation 
theory developed in \IQ in order to include the large q sector. This 
would imply to deal with the small divisors without any special Tay- 
lor expansion. Moreover, to properly cope with the cascade of modes, 
one should understand which are the actual modes excited in the q- 
breather, according to the coefficient -B^/j. 

Next, as suggested by several numerical studies, QBs appear to 
have exponential localization for any qo- One can expect that this is 
a common property of all the Lyapunov orbits in FPU chains. The 
expected differences, when changing qo, could concern the sequence 
of modes involved and, maybe, details of the exponential decay and 
stability of QBs. 

As shown above, tail resonances in the QB energy distribution ap- 
pear due to supcrharmonic almost resonances. These tail resonances 
are also observed in the evolution of natural packets, and are responsi- 
ble for the approach to equipartition, pumping energy from the packet 
core into the tail. An analytical theory would thus allow to understand 
and describe the ultimate route to equipartition, the holy grail of the 
FPU problem. 

Below we will address some of the above points. 

3 Improved perturbation theory 

In this third part we present a result based on the Poincare-Lindstedt 
(PL) perturbation scheme: for a generic go-breather, we prove the 
existence of a threshold £* , scaling with kq = j^ri and with the per- 
turbation parameter, such that for the energy density £ qo < £* the 
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averaged harmonic energies involved in the QB dynamics are exponen- 
tially localized. In the a model, for example, we get for any n > 2 

<*U < ^o« +«H)(=r^) 2 «- (1 - n)|1117 '' | ! 

7„ = C{ap n _ t f£ qo , (11) 

where g„ is the sequence of the modes excited in the QB, (Eq ) is the 
leading averaged harmonic energy 



o 

and /Z n are defined as 

Mn := — , ~P n = , max {M- 

7?„ Z=2,...,n 

This result completes the study on QB localization in Fourier space 
performed in ^H], developing a more general estimate based on the 
same PL scheme. Among the differences between the two approaches, 
we first mention the lack of an explicit local approximation for the 
smallest divisor as in Q: in fact for any n we simply factorize the 
term rj n without knowing its value. This allows us to generalize the 
previous result, once it is clear how to determine q n . Another difference 
concerns the harmonic cutoff of the leading solution ©: in order to 
perform a rigorous approach we have decided to keep all the harmonics 
in the Fourier expansion of the solution. This leads to the constant 
factor C > 1 in which takes into account both the accumulation of 
frequencies' products and all the possible non-zero higher harmonics. 
This affects the threshold £* , which is derived from (jl 1|> by requesting 
In < 1: 2 

£ *qo : = C [ — ) ' //:= max {/!„}. (13) 

We have already mentioned the scaling properties of £* , which de- 
pends on the wave number kq, but not on the actual system size. The 
main reason is that, when changing both qo and N + 1 while keep- 
ing their ratio fixed, the frequencies involved in the QB evolution are 
always the same. Then the small denominators also depend only on 
n and the same holds for /i. On the other hand, if we fix q while 
changing the size N of the system, the smallest possible divisor rjj goes 
to zero. Consequently the threshold £* g vanishes and the QB loses its 
localization in such a case. In general we are not able to provide an 
explicit value for In some cases, however, it is possible to make a 
suitable approximation similar to (J7J. For example, when go is close 
to zero or to N, it follows 

c* ~, —2 4 4 c* ~, —2 4/1 \4 
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for the a model and 

£* qo « /T Vk*, £* qa « /TV(1 - Ko ) 2 , 

for the (3 model. In these two specific cases the slope of the exponential 
decay is given by | In 72 1 : this value is preserved for almost all the 
energy distribution, until values of \ij greater than 1x2 appear in the 
tail. These higher order anomalies are the ones producing resonant tail 
peaks: the factor /i n in 1)11(1 amplifies those modes whose frequencies 
are almost resonant with uo qo . This effect is further influenced by uj nqo , 
which, depending on its value, can weaken or strengthen the amplitudes 
(see section 3.4 and Fig. |3J). An unavoidable limitation is that the 
sequence /ij is constant with respect to the perturbation parameters 
a and £ qa . Thus it does not take a possible change of position and 
intensity of the peaks in the energy distribution into account. 

In the following, we will recall the basic ingredients of the Poincare- 
Lindstedt method. Then we present the rules used to predict the 
cascade of modes excited in any (70-breather. Proposition 13.11 can be 
suitably extended to any homogeneous potential. It is strongly re- 
lated to some symmetries of the system, as explained in Appendix A 
and B (see also |32|). Finally we move to the main general statement 
about localization, performing the proof for a, j3 and a combined a — (3 
model. We finish with an application to the high-frequency-mode case 
\q a -N\<^N. 

3.1 The Poincare-Lindstedt method. 

The Poincare-Lindstedt expansion is a classical perturbation method 
used to continue a periodic orbit with respect to a small perturbation 
parameter, when fixing the amplitude (or the energy) of the system. 
The crucial idea is to expand both the solution and the modified fre- 
quency with respect to this parameter: this is essential to kill secular 
terms which appear in the recursive scheme. 

For example, consider directly the a-FPU system: we start ex- 
panding in a = . " both the normal mode coordinates and the 

V 2I\ -{-2. 

new orbit's frequency 

QS) = Xy'Q^C*). = Xy w <0. 

We define a new stretched time variable t = f2(er)i and we rewrite the 
equation of motion with respect to the new time derivative -4- 

d 2 

fL(<r) 2 —Q q (T) +uj 2 q Q q (T) = -a ^ u q u>iUJ m Qi(T)Q m (T)B q i m . 

l<l,m<N 

(14) 

If we expand (|14Jl with respect to a we get an infinite sequence of 
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equations for every term Q q 

N n-1 

UJ q ^ n B q l m YQl h) Q l m~ 1 ~ h) ^) 
l,m=l h=0 

where 

ii>0 <=0 

We need solutions of (|15|l which are 27r-periodic in r, so bounded. 
The initial condition will be chosen so to guarantee this requirement, 
and the term a;w will be used to prevent the existence of secular 
terms (unbounded motions) which could appear in the equation of the 
continued mode. 



3.2 Cascade of modes. 

The QBs computed in Fig. present an interesting feature that was 
already visible in [T2|: only a small subset of modes is taking part to 
the dynamics, the other being essentially at rest. So, for a full descrip- 
tion of such solutions, we aim at understanding which mechanism is 
responsible for this modes' selection. This is also essential in order to 
get some straightforward result from equation (|1 B|) . The main rule is 
that, at any n-step of perturbation theory, just one new mode q n is 
excited: this corresponds to the only mode in the interval [1, N] whose 
frequency equals 2 1 sin ( 2N+2 ) I • ^ ne coun t of q n stops when either 
nqo is a multiple of N + 1 or when nqo and (n — l)go have the same 
frequency (i.e. they are symmetric with respect to one of the two dis- 
persion edges). A more rigorous claim is the following (for a sketch of 
the proof see Appendix A): 

Proposition 3.1 At any k—1 perturbation order, only one new mode 
q k with frequency 2 1 sin ( 2N+2 ) I * s exc tted. This mode reads 

- f (kqo) mod(27V + 2), (kq ) mod(27V + 2) < N + 1 

qk ~ \ 2N + 2 - (kq ) mod(2iV + 2), otherwise 

(16) 

Thus it holds 

Q( m) (t) = Vm<*-1. 

The above proposition exactly predicts which modes are actually 
excited by the QB's dynamics. These orbits are strictly related to 
some of those invariant submanifolds whose existence has been proved 
in |32| using discrete symmetries of the system. However, differently 
from those arguments, the above proposition provides also an order 
according to which the modes take part in the evolution. 

The final scenario is the following. Take 1 < qo < N and go := 
gcd(qo, 2N + 2) (as in Proposition l5.1|) . For go = 1 the submanifold and 
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the QB are simply the whole mode space. For go > 2 and g ^ $ N 
the submanifold has lower dimension that the number of modes, at the 
same time it can not be embedded in a subsystem with fixed boundary 
conditions. Hence, the QB does not involve all the modes, but is not 
a rescaled solution of a lower dimensional system. For go > 2 and 
go = gcd(N + 1, go) the submanifold is embedded in a fixed boundary 
condition subsystem with -^p- particles and the QB is the ^-mode 
rescaled solution, constructed as in |25| . 

3.3 Exponential localization 

Next we prove the exponential decay of the harmonic energies for a 
generic QB. This proof makes use of the sequence of small denomi- 
nators previously denoted as T) n<m , since it naturally arises in the PL 
method. The Fourier coefficients in the PL scheme, at any new step, 
are naturally defined via a convolution product. Since I 2 is not closed 
with respect to this product, the choice of a || • > || • ||;2 norm seems 
natural to control the main physical quantities. It is then useful to 
recall the following 

Lemma 3.1 Let l\ denote the Banach space of absolutely convergent 
bi-inftnite sequences, endowed with the norm ||a|| ; i := X^jez \ a j\- Then 
l\ is a Banach algebra with respect to the convolution product 

(a*b) m = ^ a i b i 

i-\-j—m 

and 

||a*%i<||a||,i||&||^ (17) 

holds. 



3.3.1 The a model. 

We develop an estimate for the averaged harmonic energies of the lead- 
ing terms 1 ' , based on the scheme presented at the beginning of 
this section. By applying Prop. 13.11 and with the same notation for 
the modes, equation 1)15(1 reads 



dr 2 



'ik 



± 



At the first step (k = 2) 



(jfe-1) 



k-1 



E/n('-l)o(fc-l- 



1 = 1 



(18) 



cos(2t) 



' J 2q 



Ho 
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We define the quantities 

V2,l ■= - /2 ^ol' V2 := ;= min 2 {77 2i( }, M2 := -L. 
The solution becomes 

2 

= X! ^^expfmr), C 2 , m = C 2 ,- m , C 2 ,i = 0, 

ro=-2 

and the Fourier sequence satisfies 

ll&Hij < M2||C'l||fiC 2 , c 2 = UJ 2go UJ 2 qo . 

Now, we state and prove the main 

Theorem 3.1 Define for any k > 2 the following quantities: 
V k ,i ■■= \ul qa ~l 2 u 2 qQ l r] k :=^mm {q kl }, 

fi k := JI k = max { W }. 

?7fc i=2,...,fc 

Then, the leading solution Q k k qo ^ (t) reads 
k 

Q7h~ 1) ( t ) = X! C*fc,m exp (imSli), Cfc,m = 0, m = fc- 1, . . . 1 - k, 

m— — k 

and there exists a positive constant a > 1 such that the Fourier sequence 
C k fulfills the following estimate 

\\C k \\ii < Mfc (Mfc-i)^ 2 ll^illficfe, c k < w* a* -1 . 

Moreover, if £ qa — represents the energy density of the mode 

q , there exists a threshold £* a such that, for £ qa < £* o , the leading 
averaged harmonic energy fulfills 

7fc = 2 (aa/Jfc-i) £ 9o . 

proof: By induction. The case k = 2 has been already proved. Notice 
that it is not necessary to define ~p 1 , since its exponent is 0. Suppose 
the statement is true for k = n — 1 and let us prove it for k = n. 
The time-dependent force is composed of harmonics with frequencies 
running from to n 

n 

Q%~ 1 Ht)Q£-_t~ l \t)= Yl (Ci*Cn-i) m exp(imnt), 

m— — n 
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with 

(Ci * C n -{) m — (Ci * C„_;)_ m , 

and 

(Ci * C n -i) m = 0, m = n - 1, . . . , 1 - n. 

From the definition of the force Fourier coefficient 

n-l 

= X! B q n .qi^-l(Cl * Cn-l)rnUlq„U(k-l)q 
1=1 

we can expand the time dependent force 

n 

m=—n 

and write the solution as 



p(n-i) 

Qi" >(t) = -w„ 90 ^ ,, 2 ": ? I 2f , 2 exp(imOt) 



u„„ — m*u>± 

m=-n n 9a la 



= Cn iTO exp (imttt). 

m= — n 

Now we move to the estimate: 



77,1 1=1 



< 



n-l 



y n<j 



1=1 



^WlmWCn-Ofloll^llijH^n-llllJ < 



n-l 
, n-2 ' 



< -zr\\^A\\{^ n -i) n ^ui qa u: {n _ l)qa cic n -i 
Vn i=i 

= MnllCilIJi (M„-i)"~ 2 c„, 



where we have defined by recurrence 

n-l 



^^lq a ^(n-l)q ClC n -l. (19) 



;=i 

The values of the first elements 02,03,04,... suggest the possibility to 
collect a factor w™ o from any term c„: indeed, if we set 6,2 = u)2 qa and 

(n-2 
2d n _ia>( n _ 1 ) 9o + E u)i qa u)( n _i )qo did n 
J=2 
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then c n = w"d n , with d n < a n . The bound on the amplitudes 
c n < a;™ a™ , which is proved in the Appendix C, allows to estimate 
the averaged leading terms of the harmonic energies (|12[) 

2n-2 rr 1 



n 



m— — n 



is obtained by differentiating Q^™ ^ (t) with respect to time t and 
recalling that £1 = Lo qo + 0(a). The period- average yields to 



|m| <n 



|m| <n 

,1 

2 l 



< - 2 "- 2 ^« + ™ 2 <)iicyi 2 ;- 



By means of the previous estimate on ||(7„||ii we get 

VMn-1/ 

which is an exponential decay until 

ln = 2(aaJl n _ 1 ) 2 £ qo <l, Vn = 1, . . . , N - 1. (20) 

If we define (i = max/I^, the condition 

7 <1 ^ £ 9Q <£; :=2(aa M )- 2 (21) 
repesents the go specific energy's threshold for the exponential decay. 



3.3.2 The /3 model. 

For the /3-model we give the main ideas how to extend the perturbation 
estimates performed for the a model. The quartic part of the potential 
can be rewritten as 

13 N (3 

1 E ( x o+ 1 ~ x ^ = 8C/V + I) E uJiUJ-jUJhVkQiQjQhQkBijhk, 

j=0 ^ ' l<i,j,h,k<N 
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where we set B ijhk = J2 ai , 2l3 =±i ^-i+^j+^h+^k, with 

( 1, i±j±h±k = 

A i±j ± h±k = I -1, i±j±h±k = ±2(N+l) . (22) 
[ 0, otherwise 

The equation of motion for a normal mode then reads 



I) Uq^i^jUhQiQjQhBqijh- (23) 



Q q + J q Q 9 2(N + 1) — 

y ' l<i,j,h<N 

The small parameter used in the Taylor expansion is a = 2 n+2 • With 
this choice, the equation at every order s becomes 



r=0 



N 



- LO q ^ UiLOjUJhBgijh ^ Qi^QjQh- 

i.j,h—l Z+m+n— s— 1 

Following both and Proposition 13.11 we can prove 

Proposition 3.2 At any s — 1 perturbation order a new mode q s with 
frequency 2| sin 



(2s-l)<j 7T 

2N+2 



excited: 



Is = 

Thus 



(2s - l)q mod(2iV + 2), (2s - l)q mod(2iV + 2) < N + 1 
2N + 2 - (2s - l)q mod(2iV + 2), otherwise 

(24) 



'(t) = Vm<s-1. 
Consider s = 2: the solution of 



is given by 



Q«(r)=± W3?0 <A^ 



3 cos(r) 1 cos(3r) 



4 oj2„ - uj% 4 wL - 9cj2 

3?o 90 oijo "Jo 



or, after expanding cosine terms 

3 



^2 ( T ) = X! C^mexp^mfit), (72,2m = 0, C 2 , m = C 2 .- 



m— — 3 

We define the following quantities 

,2 ;2, ,2 



^2,i := K 9o ~ Z »» := min {772,;}, /i 2 := — , 
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then 



|C 2 ||,i < M2||Cx||fxC 2 , C2=W3 W W^ ) . 



Theorem 3.2 Define for any s > 2 the following quantities: 
Ms := — , M s = max 

77s l=2,...,s 

Then, the leading solution Q^ s (t) reads 

2s-l 

Q i l~ 1) ( t )= £ C Stm e X p(imQt), 4,2m =0, (25) 

m=-2s+l 

and there exists a positive constant a > 1 sitc/i that the Fourier sequence 
C s fulfills the following estimate 

ll&ll,! < Ms (M s -l)^ 2 H^illfx— 1 C„ C S < U&-V- 1 . 

Moreover, if £ qa = represents the energy density of the mode qo, 

there exists a threshold E* Q such that, for £ qa < £* o 



(E- qs )<E qo (^ 2s _ 1)qo +(2 S -lfu;l) (J±-) e-^~^ 

\Ms-l / 

with 

7s = Pafl^E^. 



proof: 

We give some hints how to modify the proof of Theorem 13.11 The 
nonlinearity in the equation of motion for QL S is 



W (2s-l) 9o J! UJ (2l-l)q UJ ( 2m -l)qo UJ (2n-l)qoQ'q l l 1 Qq 



>!' 0/^(m-l)g(n-li 

/+m+n— s+1 



where l,m,n> 1. We can define the force Fourier coefficient 



and obtain the time dependent force 

2s- 1 
fe=-2s+l 

This, as before, yields the solution 

2 S -i p^- 1 ) 

ft=-2s+l (2s-l)g 9o 
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2s-l 

= exp (ihflt). 

h=-2s+l 

The conclusion of the proof follows Theorem Kill 



3.3.3 The a - (3 model with j3 = \a 2 . 

We finally consider the a — (3 model, under the restriction = |a 2 . 
This dependence is due to the need of expanding the solution with 
respect to one perturbative parameter only. Moreover, with this choice 
the potential coincides with the fourth truncation order of the Toda 
potential and the scheme performed can be adapted to any order. 

In this case, defining a = ^ 2 n+2 ' ^ e equation at any order s — 1 
reads 



r— i)J = 1 l-\-?n—s—2 

2 ^ 

_ 3^9 E UJ i UJ j UJ hC q%0 h E Q/" Q/T • 



i,j,h—l l+m+n—s — 3 

As in the a model, the set of modes involved in the dynamics is 
selected according to Proposition 00 (see again Appendix A) 

Proposition 3.3 At any k—l perturbation order, only one new mode 

>2N+2> 



q k with frequency 2 1 sin ( 2N+2 ) I * s exc ^ e ^ ; 



3fc 



f (/cg ) mod(27V + 2), (kq ) mod(2iV + 2) < N + 1 

\27V + 2- (fcg )mod(2iV + 2), otherwise 

(26) 

Thus 

Q( m) (t) = Vm<A-l. 

Ik v ' 

The exponential decay is similar to the one obtained for the a 
model, except for the coefficient c s which is responsible for the accu- 
mulation of frequencies: 

Theorem 3.3 Define for any s > 2 the following quantities: 
Vs,l ■= \u 2 sqo - l 2 M 2 qo \, n s :=^mm {77^}, 

:= — , ~p s = max {in}. 

rj s 1=2,. ...s 

Then any solution Q^ s (t) reads 



Qls (*) = E C Sjm exp(imQt), C Sim = 0, m = s - 1, . . . , 1 - s 
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and there exists a positive constant a > 1 such that the Fourier sequence 
C s fulfills the following estimate 

\\Cshl ^ (K-iY~ 2 ll^i||?ic s , c s < u s qo a s -\ 

Moreover, if £ qo = N ^ represents the energy density of the mode qq, 
there exists a threshold £* Q such that, for £ qo < £* o , holds true 

(\><^«+^)(f) 2 ^ ,H|ln H 

\Ma-l / 

with 

7s = 2 (aa/I s _i) £ qo - 

proof: By induction we consider the step q s . The quadratic and cubic 
force are composed of terms like 

s 

Q ( V )Qi Z~ 1] = S (Ci*C m ) r exp(irnt), 

r— — s 
s 

^r 1)Q £" 1)Q £ _1) = S * ^» * cx p (* rm )' 

r— — s 

with respectively I + m = s and Z + m + n = s; observe that 

(Ci * C m ) r = (C\ * C m * C n ) r = 0, V = S - 1, S - 3, . . . , 1 - S. 

Since in both the cases the harmonics run from — s to s, it is possible 
to define the force Fourier coefficient as 

^ir^ = E 01 * C' s -0r'^ g o^( S -/) g o + 
Z=l 
s-2 

+ E * Cm * C , s _;_ m ) r o;j ( j a; mi j a;( s _;_ m ) ( j 

and expand the time dependent force 

it _1) = -".«, E Fi'-^expCirfit), 

r=— 5 

which leads to the solution 
Qi s_1) = -w S90 ^— — exp (irClt) = ^ exp {irQ.t). 

r=-s Ws «o r "<?o r= _ s 

It is remarkable that the both force contributions generate the same 
small denominators. The proof can be then easily finished, using a new 
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Figure 4: Parameters: N = 47, E = 0.1. (a): Energy of normal modes in a 
QB (open symbols) and a natural packet (filled symbols, waiting time 10 5 ) 
for qo = 47 and a = 0.25. (b): Similar to (a), but for = 0.25. (c): sequence 
Xj 1 which corresponds to case (a), (d): sequence xj 1 which corresponds to 
(b). 



definition for c s 

/ s-1 



Cs -—^sqo I ^ ^Iqo^is-^qoClCs-l + ^ UJ lqo UJ mqo^(s-l-m)q c l c mC s -l-rn 



W=l 



1=1 



3.3.4 Application: the high mode case. 

We consider the case when q$ corresponds to a large frequency, i.e. 
JV-^oCl. As for a low frequency initial excitation, also this special 
case is related to the equipartition issue. We recall, for example, that 
in |85| the energy flow from initially high-frequency-mode excitation to 
low modes was studied, restricting to a /3-PPU model. The claim was 
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that, due to the same mechanism of resonance overlapping described in 
the equipartition time T eq scales as a power law with the inverse 
of the energy density, T eq ~ £~ 2 . However, the authors didn't exclude 
also the possibility that, below a small energy density threshold, the 
dependence T eq (£) could be even slower than any power law. The 
existence of exponentially long waiting times to equipartition, scaling 
with the inverse of the energy density, and the existence of two different 
time scales for the energy spreading among modes, have been later 
numerically confirmed in |28| . 

While the impact of QBs on these observations is essentially the 
same as for the low frequency case, a new element in the structure 
of the QB itself arises. For energy densities small enough, it is still 
evident an exponential decay of harmonic energies starting from the 
initially excited mode. However, a remarkable difference with the low- 
mode case exists and it is related to the effect of the a term in the 
potential: if in the pure model the energy is mainly shared among 
high frequency modes, in the a case the energy distribution involves, 
from the very beginning of the evolution, also low frequency modes. We 
compute QBs with high frequencies, and plot their energy distributions 
in Fig.0] The panels (a) and (b) of Fig. 0] show two g-breathers with 
the same largest q , for the a and (3 model, respectively. We also 
plot the energy distribution of the corresponding natural packets. We 
again recover that both distributions are remarkably similar. Note that 
for the /3-model the modes decay monotonously with decreasing mode 
number. However for the a-model a strong excitation of low frequency 
modes is observed as well. 

The two lower panels (c) and (d) in Fig. ^represent the Xj sequences 
corresponding to the two QB periodic orbits. Since in the forthcoming 
part we want to obtain, from the general definition of £* in Theorem 
13.11 an approximated but explicit dependence of the threshold energy 
density on kq := -^j, we will actually consider the sequence r\j in- 
stead, since rjj is still very close to Xj- We deal with the two initial 
small divisors 772 and 773, arising respectively from the first and second 
PL step. In the first step, the low mode 2N + 2 — 2q is coupled to 
qo, but the forcing term Q q ° ' (t) includes only even harmonics, thus 
generating in the solution a small denominator — ^— . We expand u>2q a 

around 7r and obtain the approximation 

"l qo = 4 sin 2 (tt(1 - no)) « 4^ 2 (1 - k ) 2 . 

In the second step, instead, the mode 3go — 27V — 2 is excited. The 
odd harmonics of the time periodic force generate the small denom- 
inator — _ 2 ■ Since both uj% qo and uj qo are very close to 2, their 

difference is expected to be quadratic in (1 — kq). We expand u>3 qa in 
a neighbourhood of §7r 

3 3 3 9 

uj 3qo = 2 sin (-7tk ) = 2 cos (-7T - -7tk ) = 2 - ttt 2 (1 - k ) 2 + h.o.t., 
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and oj qo around ^tt 

u qo = 2sin(i7TK ) = 2cos(i7r - ^vtko) = 2 - ^tt 2 (1 - «o) 2 + h.o.t., 
thus deducing 

K D -<l«87r 2 (l-K ) 2 . (27) 

So, in agreement with panel (c) in Fig. we have finally explained 
the relationship 7/2 > Vs- Thus we may expect (Eq 2 ) > (Eq 3 )- To 
compare panel (c) with panel (a) in Fig. 01 we recall that the averaged 
harmonic energies (Eq B ) are also corrected by the coefficients (wJ 9o + 
s 2 ujq ) ■ Finally, from the previous local approximation of strong nearby 
resonances, one can deduce 

£;«a~V(l- Ko ) 4 . 

It is straightforward to deal with the (3 case and obtain from Theorem 

£;«/rV(i- Ko ) 2 . 
4 Conclusions. 

We have studied the energy localization of a generic Lyapunov orbit, 
called g-breather, in the pure a and /3 FPU model. We have proved, 
thus generalizing a result from |16| . that any ^-breather exhibits an 
exponential localization in the Fourier space if the qo-mode energy den- 
sity £ qo is below a threshold £* Q which depends only on the coupling 
parameter (either a or /?) and on the wave number kq = ^° 1 . The 
specific value of £* Q suggests that a /3-QB will delocalize for higher 
energies. Moreover, the exponential decay of amplitudes in a QB runs 
through a subset of normal modes: we revealed the algebraic rule which 
selects the modes involved in the QB dynamics, thus explaining their 
order in the perturbation decay. This rule is strictly linked to the dis- 
crete symmetries of the system. Any QB lies on one of those invariant 
submanifolds described, for example, in and 

The possibility to evaluate the Poincare-Lindstedt (PL) perturba- 
tion scheme at any order, allows a qualitative explanation of the previ- 
ously observed resonant peaks in the tail of the energy distribution. In 
the numerical experiments on natural packets with initial low frequency 
modes, we notice anomalies in the tail of the energy distribution. They 
pump energy from the natural packet core into the tail. We think that 
this effect is of central relevance to the pathway to equipartition. Such 
anomalies are also characteristic of any QB, even if in the low fequency 
mode simulations (see Fig.n an d Fig. 01 they are more visible. While 
these peaks have been plotted (but not discussed) already in JS], we 
have developed a scheme which is capable of explaining the positions of 
these peaks. We have shown that ther is a direct connection between 
the position and intensity of these peaks and nonlinear resonances in- 
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volving the QB's frequency Q, 



X„,m=Ko- m2 °«ol- (28) 

In the limit of very small perturbations, the sequence (|28|) reduces to 
Vn,m = \u 2 nqo ~ m 2 Uq \, (29) 

which is the one arising in the PL scheme. We have used the latter 
one for a rigorous analytical approach. We still lack a suitable Fourier 
expansion which allows to implement recursive estimates including the 
true small divisors (|28fl. Nevertheless it is possible, restricting to low 
and high frequencies, to improve the estimate of Theorems 13.11 and 
13.21 with a modified PL approach, thus discovering that f3 QB thresh- 
olds are larger than a ones. Finally, we have also extended the one- 
parameter PL perturbation scheme to a special a — (3 model, with 
(3 = | a 2 . Thus the PL method developed, although fitting the numer- 
ical simulations only for perturbations small enough, gives a rigorous 
analytical way to describe localization in Fourier space for any homo- 
geneous FPU-like potential and any truncated Toda potential. 

We did not consider the linear stability fo QBs here. The reason is 
that we focussed on the exponential localization of QBs. This property 
we believe is not depending on whether the QB is stable or not. Let us 
explain this point. If a QB is exponentially localized, the core can not 
be in strong resonance with the tail modes - otherwise the QB would 
not localize. Thus almost any instability which a QB may encounter 
on the path of its continuation, will affect the dynamics of the modes 
in the core of the excitation. Indeed, it has been observed in ^U], that 
such an instability triggers weak chaos among the core modes, leaving 
the localization intact. The threshold to this weak chaos state was 
correctly identified as an instability threshold of QBs Thus, given 
a QB is localized, its stability or instability will not change the fact 
that nearby trajectories (natural packets) will stay for long times in 
that phase space part. 
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5 Appendix. 



5.1 Appendix A. 

We give here a sketch of the proof of Prop. 13.11 We have decided to 
omit the proofs of Prop. 13.21 and Prop. \'6.'6\ since, a part from small 
variations, they are based on the same ideas. 



Proof of Proposition I3.lt We give an inductive proof. For 
k — 2 the statement is true. Suppose, by inductive hypothesis, that 
it is also true until k — n — 1. Let us introduce here the following 
notation for the modes excited qo — q 1 ,q 2 , ■ ■ ■ >? n _i an d the repeated 
dispersion law F(x) :— 2| sin( 2 j^_ 2 )|; the frequencies are respectively 

F(q ), F(2q ), . . . , F((n - l)q Q ). 

By the inductive hypothesis, any mode q~j is obtained from jq$ through 

f J(2N + 2)-jq , F'(jq o )>0, 
q i \jq -J(2N + 2), F'(jq o )<0, 

for some J related to jqo. From 1)15(1 we know that the equations of 
motion at the next order include products like 

r__ n(°W™^ 2 ) 

B-- Q (1) G> ( "~ 3) 



and, among them, the only ones giving a new frequency are 



)(0) n (n-2) 

->(lWn-3) 

'92 ^9„-2 



Indeed, by induction Q^ n h x ' = for m > n — h, and q ± q t ± q m = 
mod(2iV + 2) is a necessary condition for B q q i q m ^ 0. We claim that 
the equations giving the new q modes 

q±q 1 ±q n _ 1 = 0,2N + 2, 
q±q 2 ±q n _ 2 = 0,2N + 2, 
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all admit a common solution, which reads 

= (N(2N + 2)-nq , F'(nq )>0, 
qn \nq a -N(2N + 2), F'(nq )<0. 

As an example we consider just the second equation in the previous 
list. Let us assume that 

q 2 = l(2N + 2)-2q , 
q n _ 2 = L{2N + 2)-(n-2)q , 

then the possible combinations are 

?2 + 9«-2 = (L + l)(2N + 2)-nq , 

-Q2~Qn-2 = -(L + l)(2N + 2) + nq , 

~h-~la-2 = (I- L)(2N + 2) + nq -4q Q , 

-?2 + 9«-2 = (L-l){2N + 2)-nq + 4q Q . 

The last couple give the already included mode g„_ 4 , so we focus on 
the first couple, giving the four equations 

Q + q 2 + Qn-2 = 0, q - q 2 - q n _ 2 = 0, 

q + q 2 + Qn-2 = 27V + 2, q - q 2 - q n _ 2 = 2N + 2; 

notice that the first and the last above equations don't admit a solution 
q G [1, N], so the only solutions are 

q = 2N + 2-q 2 -q n _ 2 = nq -(L + l-l)(2N + 2), 
Q = q 2 + q„-2 = (L + l)(2N + 2)-nq , 

when respectively q 2 + q n - 2 > N + l and q 2 +q n _ 2 < N + 1 ; they both 
correspond to the same frequency 2 \ sin ( 2N+2 ) I • 

■ 

5.2 Appendix B. 

Let 

R : (Q j ,P j )~(Q j+1 ,P j+1 ) 

S : (Qj, Pj) 1— > (— Q2N+2-3, — P2N+2-]) 

the discrete symmetries of the FPU system and, for a generic symmetry 
G on a manifold M, 

Fix(G) := {x G M I G(x) = x}. 
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It is possible to prove the following 



Proposition 5.1 Let qo be a mode to be continued, go := gcd(2N + 
2, go) & n d K o '■= 2N+2 . Then, the qo-breather sequence of modes 
coincides with Fix(R K °) n Fix(S) considered as an invariant subman- 
ifold of Fix(S). 

proof: Firstly, following 32 , we recall that 

Fix(R K ° ) = {Q 3 ; = Pj = iffjV Omod(go)}, 
Fix(S) = {Qj = -Q 2N+2 -j, Pj = -P2N+2-j}, 

so that it is possible to interpret Fix(S) as the embedding of the fixed 
boundary condition system in the periodic boundary condition system. 
The thesis reduces to 

{j = Omod( 5o )} = {q k }, j = l,...,N. 

In order to show the double inclusion of the two sets it is useful to 
remind that, from the definition of go, it follows 

aq +f3(2N + 2) = ±1 a,(3eZ. (30) 

C) let us take j = 0mod(<7o)j which means j = Jgo with J < [~^]; 
we have that 

Aq + B{2N + 2) =±Jg , A — aJ, B — [3J. 
2) this case is trivial, since the definition in Proposition 13. II implies 

milk- 



Proposition ^ . II claims that we should interpret the fixed boundary 
conditions system as embedded in the periodic one and then consider 
the symmetries Fix(R K ) with 2(7V + 1) = Kg; observe that it must be 
k^2, otherwise go = N + 1 > N. So, any QB is a special orbit of one 
of such submanifolds. Observe also that, for k = 2N + 2 the invariant 
submanifold coincides with the fixed boundary subsystem Fix(R K ) D 
Fix(S) = Fix(S), while for k = 2 we have Fix(R K ) n Fix(S) = {0}. 
All the others submanifolds Fix(R K ) n Fix(S) with k not being a 
divisor of 2N + 2 are meaningless since 

Proposition 5.2 For all I = 1, ... , 2N + 2 holds true 

Fix{R l ) = Fix(R scd{L2N+2] ). (31) 
proof: Simply, by calling re = gcd(Z, 2N + 2) and 2N + 2 = ng 

2N + 2 _ 2N + 2 

gcd(Z,2iV + 2) ~ 9 ~ gcd(K, 2N + 2) ' 
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The final scenario is the following. Take 1 < go < N and go, kq ad 
defined in the Proposition l5.il then 

1. g Q = 1 ^ k = 2N + 2 O Fix(R K «) n Fia;(5) = Fix(S): the 
submanifold is trivial. 

2. on > 2 and 2N+2 d N =>. the submanifold doesn't coincide with 
Fix(S) but cannot represent a fixed boundary conditions subsys- 
tem. Hence, the QB does not involve all the modes (not trivial) , 
but it is not a rescaled solution of a lower dimensional system. 

3- go > 2 and go = gcd(N + 1, on) =^ the submanifold is the embed- 
ding of a fixed boundary condition subsystem with N ^ particles 
and the QB is the —-mode rescaled solution, constructed as in 

m 

5.3 Appendix C. 

It is known (cfr. [2]]) that the sequence 



Ai = 1, 

A„ = Ysj=l AjA„_j = Ysi+j=n AiAj 



(32) 



is equivalent to a suitable functional equation 

g(z) = z + g 2 (z), g{z)^Y. X - zr ( 33 ) 



and that follows X r < 4 r 1 . Consider now the modified problem 

Ai = Hi, 

A n = /i T i Zjj=i AjA„_j ' 



(34) 



where fj,j is a limited sequence. Let /i := sup^ /ij and define the auxil- 
iary sequence 

whose related function satisfies 

g(z) = fiz + ng 2 (z). (36) 
It is easy to prove by induction that 

n— 1 n— 1 

A n = /i n ^ ' Aj An-j < /i ^ ' AjA„_j = A„. 
3=1 3=1 

Thus, a bound for A„ can be obtained from the following asymptotic 
behaviour of A n 



2^-1 ( 2 ™_3)!! M^.,2 W -i 



A„ = /i(2/i 2 )™ -1 ^— — -j— ^— — -(4^ 2 ) 
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This result can be applied to Theorem 13.11 where we deal with the 
sequence 

C n := UJ nqa ljJ i<;o ljJ (n-i)9o CiC ™- ( ' 
2 = 1 

which reads also 

n-i 

An ■= Un qo hK-l, 

2=1 

where A; = LUi qo ci and u>f < 4. 

Let us move to a more complex sequence 



Ai = 1, 

An = Si+i+fe=n+l AjAjAfc 



(37) 



where plainly 1 < k < n — 1. Following 2 , we define an ana- 
lytic function g(z) — X)r>i Ar^ r such that (|37|l is equivalent to the 
functional equation 

g = z+- z g 3 ; (38) 

we stress that, since g has a first order zero in z = 0, then also 
\g{z) is analytic. Let us define an auxiliary function f{z) = ~g(z) = 
^ r>0 c r z r , with c r = A,+i, and rewrite l|38(l as 

f = l + zf, /(0) = c = Ai = 1. (39) 

The differentiation of Q39J1 gives 

f 3 

1 ~ 1 - 3z/ 2 ' 

which, exploiting again (|39|) . reads as a z-independent expression 

^ = 3^7- (40) 

We may assume the following general description for J^" 1 ) 

: (3 - 2/)i ' 



2r-3 4r f fc 



/<-.) = ££ ft /_. (41) 



j=0 fc=0 



with suitable values for and check its reliability by induction. The 
differentiation of the generic addendum in (|41|) gives 



f k kf k+3 2jf k + 4 



dz{3-2f)j (3-2/)i+! (3 - 2f)i+ 2 ' 

which yields 

2r-14r+4 .j. 

f (r) = E E T^JT (42) 
j=0 fc=0 v J ' 
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where the new coefficients otjk are properly defined through /3jk- As 
noticed in [2], Lemma 11.2, we can avoid the computations of ctjk 
simply using the previous computation 



2r-3 4r 



fWrm VV/i fc/(0) i 2jm 



_ it (3-2/(0))^i (3-2/(0))^ 



2r-3 4r 



j=0 fe=0 



which allows to find 

2r — 3 4r 

C ^ = < ^ E Eft* = 8ri/ (r_1) (0) = Scr-i- 

j=0 fc=0 

This concludes the proof, since we have finally shown 

Let us consider the sequence 
Mi = 1, 

m = 1, . (43) 

Here the solution comes from the comparison of (|43H with the previous 
sequence IpTfjl. Indeed, with the choice Ai = /Xi = = 1, it is possible 
to prove by induction that ji n < X n -i < A„. 
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List of captions 



1. FIG.l: (a): Natural packet evolution initially exciting mode 
q = 1. The energies of normal modes are plotted versus q. 
Circles - 10 4 , squares - 10 5 , rhombs - 10 6 . Dashed line - g-brcather 
from (b) for comparison, (b): The energies of normal modes 
versus q for a q-breather with q = 1. (c): sequence r/J 1 . (d): 
sequence xj 1 - Parameters in all cases: N = 31, a = 0.33, £ = 
0.01. 

2. FIG. 2: (a): Natural packet evolution initially exciting mode 
qo — 2 (circles) . The energies of normal modes are plotted versus 
q at t = 10 5 . The dashed line shows the energy distribution for 
a corresponding g-breather with qo = 2. (b): Same as (a), but 
for (jo = 3. (c): sequence rjj 1 (2,47). (d): sequence r)J 1 (3,47). 
Parameters in all cases: N = 47, a = 0.1, £ = 0.0025. 

3. FIG.3: Parameters: N = 47, a = 0.1, £ = 0.0025. (a): Energy 
of normal modes versus mode number of a QB for q = 2. (b): 
Similar to (a), but for a QB with q = 3. (c): sequence xj 1 (2,47). 
(d): sequence xj 1 (3, 47). 

4. FIG. 4: Parameters: N — 47, £ = 0.1. (a): Energy of normal 
modes in a QB (open symbols) and a natural packet (filled sym- 
bols, waiting time 10 5 ) for qo = 47 and a = 0.25. (b): Similar 
to (a), but for (3 = 0.25. (c): sequence xj 1 which corresponds to 
case (a), (d): sequence xj 1 which corresponds to (b). 



38 



